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Abstract.  Various  features  of  solar-sector  synchronous  modulations  of  the 
particulate  cosmic  radiation  reaching  the  earth's  atmosphere  have  been  stu¬ 
died  using  satellite  and  surface  data.  The  flux  in  the  broad  maximum  of  the 
galactic  cosmic  ray  differential  spectrum  (near  1  Gv  rigidity)  exhibits  an 
intermittent  north-south  asymmetry  in  mid  and  high  geomagnetic  latitudes. 
This  modulation  exhibited  a  strong  association  with  geomagnetic  disturbance 
index  and  interplanetary  magnetic  field  direction  during  the  1 96 4  and  1965 
years  of  sunspot  minimum.  Such  correlations  are  consistent  with  the  predic¬ 
tions  of  a  theory  that  attributes  the  north-south  asymmetry  to  reconnection 
of  the  interplanetary  and  geomagnetic  fields.  This  finding  is  also  con¬ 
sistent  with  suggestions  that  solar  activity  influence  on  atmospheric 
processes  may  be  mediated  by  the  resulting  modulations  of  upper  tropospheric 
ionization. 

Introduction 

Two  modulations  of  galactic  cosmic  rays  at  earth,  discovered  in  1967 
satellite  data,  have  been  previously  described  according  to  a  magnetic  cou¬ 
pling  model.  That  model  explained  the  modulations  as  due  to  combined  effects 
of  the  interplanetary  and  geomagnetic  fields  (IMF  and  GMF)  (Ely,  1977).  One 
of  the  modulations  was  a  north-south  asymmetry  (USA)  in  which  the  flux  from 
400  Mev  (rigidity  ca.  1  Gv)  to  over  1  Gev/nucleon  was  observed  to  decrease 
approximately  30 %  in  mid  to  high  northern  geomagnetic  latitudes  during  posi¬ 
tive  sectors  of  the  interplanetary  magnetic  field  (IMF  away  from  the  sun), 
and  in  southern  latitudes  during  negative  sectors  (IMF  towards  the  sun).  The 
NSA  was  ascribed  to  reconnection  between  the  GMF  and  the  IMF,  in  that  the 
galactic  cosmic  ray  flux  is  expected  to  decrease  rapidly  in  the  flux  tube 
between  the  connected  pole  and  the  sun  (since  the  sun  is  not  a  significant 
source  at  these  energies). 

In  the  other  modulation,  the  equatorial  flux  also  exhibited  an  IMF  sector 
polarity  dependent  variation  of  similar  amplitude,  but  showed  a  seasonal 
character  with  opposite  phases  in  summer  and  winter.  That  modulation  was 
explained  to  result  from  parallel  and  antiparallel  configurations  of  the  IMF 
and  GMF  affecting  cutoff  rigidities.  Both  modulations  exhibited  similar 
rigidity  dependence  and  were  found  at  about  IX  amplitude  ca.  10  Gv  in  cosmic 
ray  neutron  monitor  (NM)  data.  The  smaller  amplitudes  seen  in  the  NM  records 
are  consistent  with  the  higher  mean  rigidity  of  primaries  contributing  to 
those  data.  Both  modulations  appeared  to  exhibit  complex  behavior,  the  NSA 
amplitude  becoming  maximum  and  the  equatorial  modulation  disappearing  com¬ 
pletely  near  the  1964-1965  sunspot  minimum.  The  NSA  appeared  to  be  associ¬ 
ated  with  increases  in  the  geomagnetic  disturbance  index,  ap,  consistent  with 
the  GMF-IMF  reconnection  etiology  proposed  by  the  magnetic  coupling  model. 
To  more  precisely  define  the  association  of  NSA  with  ap  elevation  is  the 
principal  purpose  of  the  work  reported  here.  In  essence,  the  magnetic  cou¬ 
pling  model  suggests:  (1)  the  cosmic  ray  flux  from  mid  to  high  latitudes  (in 
regions  including  the  GMF  cusps)  should  be  maximum  when  GMF-IMF  reconnection 
is  not  taking  place  and  the  GMF  is  undisturbed  (i.e.  ap  near  0);  and  (2)  over 
some  range  of  ap,  the  cosmic  ray  flux  will  decrease  near  latitudes  connected 
to  the  sun  (northern  in  positive  and  southern  in  negative  IMF  sectors)  when 
reconnection  .is.  occurring.  Testing  of  this  two  part  hypothesis  is  the  major 


formal  goal  of  the  work  reported  here.  The  magnetic  coupling  model  does  not 
predict  the  more  complex  rigidity  dependent  effects  at  the  unconnected  poles. 
The  small  amount  of  satellite  data  evaluated  near  1  Gv  did  not  show  signifi¬ 
cant  change  in  cosmic  ray  counting  rate  in  the  unconnected  region.  The  some¬ 
what  higher  rigidity  neutron  monitor  data  examined  here  are  seen  to  exhibit  a 
distinct  but  significantly  smaller  decrement  at  the  unconnected  pole. 
Analysis  of  this  complex  issue  is  not  the  goal  of  the  present  paper,  and  com¬ 
ment  is  merely  incidental. 

Interest  in  these  modulations  and  in  the  resulting  changes  in  upper  tropo¬ 
sphere  ionization  is  largely  motivated  by  their  possible  effects  on  atmos¬ 
pheric  processes.  The  existence  of  such  changes  over  the  eleven-year  solar 
cycle  and  their  possible  importance  to  the  sun-weather  problem  were  pointed 
out  by  Ney  (1959).  The  existence  also  of  short  term  changes  (on  a  scale  of 
days)  as  well  as  mechanisms  to  explain  the  influence  of  all  such  ionization 
changes  on  climate,  weather,  high  cirrus  cloud  formation,  droughts,  cyclonic 
vorticity  and  lightning  frequency  have  also  been  described  (Ely,  1977,  1 97 9 » 
1984). 

Correlations  have  been  reported  between  most  of  these  and  solar  activity 
with  high  (>3  sigma)  statistical  significance  (Eddy,  1976,  1979;  Stuiver, 
1961;  Mitchell  et  al,  1977;  Roberts  and  Olson,  1973;  Wilcox  et  al,  1974; 
Hines  and  Halevy,  1977;  Lethbridge,  1981;  Scherrer,  1979;  Herman  and  Gold¬ 
berg,  1978).  Because  of  interest  in  the  possible  role  of  the  modulations 
studied  here  as  mediators  of  the  solar  activity  influence  on  atmospheric 
processes,  several  analyses  are  in  progress  utilizing  more  data  than  were 
available  for  the  earlier  study.  The  importance  of  tropospheric  ionization 
by  cosmic  rays  as  a  possible  mediator  of  the  sun-weather  correlations  is 
amplified  by  the  paucity  of  possible  candidates  to  perform  the  coupling 
(Dessler,  1975). 

This  paper  presents  an  analysis  of  the  north-south  asymmetry  modulation 
during  the  1964-1965  two  year  period  including  sunspot  minimum.  For  this 
interval,  we  show  below  that  the  NSA  exists  and  exhibits  a  highly  significant 
association  (p<<.001)  with  the  magnetic  disturbance,  index,  ap,  as  the  recon¬ 
nection  hypothesis  suggests. 

Data  Selection  and  Analysis 

The  NSA  modulation  was  first  observed  in  1967  satellite  data  and  then  con¬ 
firmed  in  NM  data  recorded  from. 1964  through  1968.  The  NSA  was  found  to  have 
a  large  amplitude  in  the  sunspot  minimum  years  1964  and  1965.  The  satellite 
(0V1-86)  measured  cosmic  rays  above  400  Mev/nucleon,  and  provided  information 
on  the  radiation  producing  most  of  the  ionization  in  the  upper  troposphere. 
Although  a  great  deal  of  data  exists  for  lower  energy  proton  and  alpha  fluxes 
(up  to  a  few  Mev),  very  little  have  been  recorded  at  the  energies  above  400 
Mev/nucleon  (with  ranges  above  200  grams/cm2  necessary  to  reach  the  upper 
troposphere)  on  low  altitude  polar  orbits  (necessary  to  cover  the  earth  and 
define  the  atmospheric  region  "illuminated"  through  the  instrument  aperture). 
Surface  neutron  monitors  (NM)  provide  continuous  records  of  cosmic  ray  fluxes 
in  the  geomagnetic  latitude  regions  of  maximum  NSA  indicated  by  the  satellite 
data.  Because  high  energy  primaries  dominate  the  NM  records  (via  atmospheric 
shielding  and  collision  multiplicity  effects),  approximately  half  the  counts 


are  due  to  energies  above  10  Gev,  and  the  NSA  amplitude  is  usually  less  than 
1$  in  these  data.  Nevertheless,  the  large  quantity  of  continuous  NM  records 
provide  a  means  of  readily  examining  the  ap  association  of  NSA  near  the 
1964-1965  sunspot  minimum,  a  time  when  the  other  cited  evidence  suggests  the 
NSA  was  present.  The  additional  studies  in  progress  are  designed  to  define 
the  NSA  behavior  up  to  the  present  time. 

For  these  reasons,  northern  (Deep  River)  and  southern  (Mt.  Wellington)  NM 
records  for  the  2  year  period  from  January  1,  1964  through  December  31.  1965 
as  well  as  ap  data  were  obtained  on  magnetic  tape  from  the  NOAA  World  Data 
Center-A  at  Boulder.  The  IMF  sector  polarities  for  the  same  years  were 
obtained  from  the  Stanford  University  Center  for  Space  Science  and  Astrophy¬ 
sics.  Data  processing  was  accomplished  on  a  VAX  11-780.  The  input  data  were 
read  from  tape  to  disc  where  NM  records  for  both  stations  were  binned  as  nor¬ 
malized  hourly  count  rates  according  to  IMF  sector  polarity  (positive  or 
negative)  and  according  to  ap  amplitude.  Normalization  was  accomplished  for 
each  hour  of  NM  data  through  division  by  the  most  recent  hourly  count  rate  at 
which  ap  was  near  zero  (to  remove  secular  trends  such  as  sunspot  cycle 
changes).  Of  the  17,161  hours  of  data  processed  for  each  NM,  86%  occurred  at 
values  of  ap  <  15  and  are  included  in  the  figures.  For  ap  <  2  no  significant 
NSA  was  detected.  From  ap=2  to  15 »  the  NM  decrement  (modulation)  increased 
approximately  linearly  and  exhibited  the  NSA.  For  the  remaining  1 4%  of  the 
data,  ap  ranged  from  15  to  179.  These  highly  disturbed  data  have  not  been 
included  here  because  of  very  poor  statistics  and  the  probable  dominance  of 
other  effects  than  "quasi-linear"  reconnection.  The  resulting  representa¬ 
tions  of  the  raw  data  are  displayed  as  Figures  1  and  2.  To  produce  Figures  3 
and  4,  the  ordinates  of  Figures  1  and  2  were  filtered  by  calculating  weighted 
moving  averages  of  order  3  (Spiegel,  1961)  as  follows:  one  half  the  NM  rate 
(ordinate)  at  each  value  of  ap  was  added  to  one  fourth  the  NM  rates  for  both 
the  previous  and  following  values  of  ap. 


Results  and  Discussion 


The  magnetic  coupling  model  predicts  that  the  galactic  cosmic  ray  counting 
rate  will  decrease  (and,  conversely,  that  solar  cosmic  rays  will  increase)  at 
the  top  of  the  atmosphere  in  regions  connected  to  the  sun  by  IMF-GMF  merging. 
In  this  and  related  papers,  we  define  the  amplitude  of  the  resulting  NSA  for 
a  given  NM  station  (at  each  ap  value)  as  the  difference  between  the  counting 
rates  in  positive  and  negative  IMF  sectors  (NM(+)-NM(-) ).  Each  of  the  four 
figures  represents  over  17»000  hours  (2  years)  of  data  recorded  at  one  of  two 
NM  stations.  The  stations  differ  in  two  important  respects.  At  Mt.  Welling¬ 
ton,  the  count  rate  is  30  times  lower  (6  x  104  per  hour  vs.  1.8  x  106  at 
Deep  River),  and  the  cut-off  rigidity  is  nearly  2  times  higher  (1.89  Gv  vs. 
1.02).  Both  of  these  differences  reduce  the  NSA  amplitude  at  Mt.  Wellington. 
The  maximum  standard  deviation  (sigma)  for  counting  statistics  in  the  ap 
range  of  interest  (3  to  15)  at  Deep  River  is  .00004,  five  times  smaller  than 
the  corresponding  .0002  value  at  Mt.  Wellington.  As  a  result,  the  separation 
between  the  positive  and  negative  IMF  sector  data  is  over  25  sigma  for  the 
former  station  but  ranges  from  only  2  to  8  sigma  for  the  same  8  ap  values  at 
the  latter.  Nevertheless,  even  at  Mt.  Wellington  the  separation  (i.e.  the 
NSA)  is  highly  significant  (p  <<  .001)  from  ap  =  3  to  15  inclusive,  as  can  be 
seen  by  forming  the  product  of  8  probabilities,  the  smallest  of  which  is  p  < 
.05.  We  have  not  yet  determined  the  cause  of  the  anomalous  transient  in  the 
negative  IMF  sector  data  at  Deep  River  for  ap  =  6  and  7.  However,  its 
existence  does  not  affect  our  conclusions. 

In  addition  to  the  statistical  noise  due  to  low  counting  rate,  the  NSA  was 
also  reduced  by  another  effect  at  Mt.  Wellington,  as  can  be  seen  in  the  fil¬ 
tered  data  of  Figure  4.  We  suggest  this  is  due  to  the  higher  cut-off  rigi¬ 
dity  because  the  separation  (NSA)  must  approach  zero  as  rigidity  increases 
without  limit.  The  rigidity  dependence  of  the  NSA  amplitude  has  been  summar¬ 
ized  elsewhere  (Ely,  1977);  its  steeply  falling  character  is  seen  from  the 
30J  value  at  1  Gv  (satellite)  and  the  0.1J  seen  in  the  2  to  10  Gv  range  (Mt. 
Wellington  NM). 

In  the  satellite  data,  there  was  no  significant  decrease  observed  in  the 
count  ra^e  for  the  unconnected  regions  (north  in  IMF  negative,  and  south  in 
IMF  positive  sectors).  This  "complete"  asymmetry  at  1  Gv  facilitates  the 
magnetic  coupling  model  explanations  of  several  atmospheric  processes  includ¬ 
ing  the  correlations  of  lightning  incidence  and  of  cyclonic  vorticity  with 
respect  to  IMF  sector  boundary  passage.  In  the  NM  records,  however,  a 
decrease  in  count  rate  is  seen  in  both  the  connected  and  unconnected  data 
(differing  in  amplitude,  of  course,  by  the  NSA).  As  mentioned  earlier,  dis¬ 
cussion  of  the  decreased  asymmetry  in  the  higher  rigidity  NM  data  (due  to 
effects  at  the  unconnected  pole)  is  beyond  the  scope  of  this  paper  but  has  no 
bearing  on  the  conclusions  presented  here. 

Conclusions 

In  the  data  sets  analyzed  here,  the  NSA:  (1)  is  not  significant  in  raw 
data  during  quiet  times  (ap  <  2;  17J  of  the  data);  and  (2)  is  present  for  ap 
values  from  3  to  15,  inclusive,  as  a  highly  significant  difference  in  the  IMF 
positive  and  IMF  negative  count  rates  (p  <<  .001).  This  work  constitutes  a 
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positive  test  of  the  magnetic  coupling  model  predictions  and  hence  is  suppor¬ 
tive  of  a  theory  which  offers  explanations  for  coupling  between  solar 
activity  and  atmospheric  processes. 

Acknowledgements .  We  thank  Drs.  J.  J.  Lord,  G.  K.  Parks,  and  R.  J.  Wilkes 
for  helpful  comments.  The  generous  sharing  of  data  and  opinions  by  Drs.  J. 
M.  Wilcox  (prior  to  1983),  P.  H.  Schemer,  and  J.  T.  Hoeksema  of  the  Stan¬ 
ford  University  Center  for  Space  Science  and  Astrophysics  are  appreciated. 
The  skilled  assistance  of  Vi  Miller  of  NOAA's  World  Data  Center-A  (Boulder) 
has  long  been  of  great  value.  The  work  was  supported  by  USRA  Contract  NAS8- 
35530,  and  ONR  Contract  N00014-77-C-0392. 


References 


Dessler,  A.  J. ,  Some  problems  in  coupling  solar  activity  to  meteorological 
phenomena,  NASA  Special  Pub.  366.  187-197,  1975. 

Eddy,  J.  A.,  The  Maunder  minimum,  Science.  192.  1189-1202,  1976. 

Eddy,  J.  A.,  Climate  and  the  changing  sun,  Climate  Change.  1.  173-190,  1977. 

Ely,  J.  T.  A. ,  Equatorial  modulation  and  north-south  asymmetry  of  galactic 
cosmic  rays  due  to  the  interplanetary  magnetic  field, 
jL-Gsoptiva.  Res.,.,..  82.  3643-3648,  1977. 

Ely,  J.  T.  A.,  Solar  activity  and  terrestrial  weather:  the  magnetic  coupling 
model,  NASA  Conference  Pub.  2098.  A  25-26,  1979. 

Ely,  J.  T.  A.,  Cosmic  rays,  solar  activity,  magnetic  coupling,  and  lightning 
incidence,  NASA  Contractor  Report  3812.  1- 1-4 ,  1 984 . 

Herman,  J.  R. ,  and  R.  A.  Goldberg,  Sun.  Weather,  and  Climate.  MSA  SP-426 . 
1978. 

Hines,  C.  0.,  and  I.  Halevy,  On  the  reality  and  nature  of  a  certain  sun- 
weather  correlation,  J.  Atmos.  Sci. .  R4 .  382-404,  1977. 

Lethbridge,  M.  D. ,  Cosmic  rays  and  thunderstorm  frequency, 
Geophvs.  Res.  Lett..  8.  521-522,  1981. 

Mitchell,  1.  M. ,  C.  W.  Stockton,  and  D.  M.  Meko,  Drought  cycles  in  the  United 
States  and  their  relation  to  sunspot  cycle  since  1700  A.D.,  EOS.  58.  694, 
1977. 

Ney,  E.  P. ,  Cosmic  radiation  and  the  weather,  Nature.  183.  451-452,  1959- 

Roberts,  W.  0.,  and  R.  H.  Olson,  Geomagnetic  storms  and  winter-time  300-mb 
trough  development  in  the  North  Pacific-North  America  area, 
J,  A taos.  Sci..  3,0.  135-140,  197  3- 

Schemer,  P.  H. ,  Solar  variability  and  terrestrial  weather, 
Rev.  Geophvs.  Space  Phvs..  17.  724-731.  1979- 

Spiegel,  M.  R. ,  Theory  and  Problems  of  Statistics  (Schaum's  Outline  Series), 
359  PP.,  McGraw-Hill,  New  York,  1961. 

Stuiver,  M. ,  Variations  in  radiocarbon  concentration  and  sunspot  activity,  j l. 
Geophvs.  »  ££»  273-276,  1961. 

Wilcox,  J.  M. ,  P.  H.  Schemer,  L.  Svalgaard,  W.  0.  Roberts,  R.  H.  Olson,  and 
R.  L.  Jenne,  Influence  of  solar  magnetic  sector  structure  on  terrestrial 
atmospheric  vorticity,  J..  Atmos.  Sci. .  31,  581-588,  1974. 


-5- 


Wellington 


♦mf  -  o 

-nrr  ■  « 

•  <  0.00004 


■  <  0.0003 


Figure  1.  Deep  River  neutron  aonltor  (WO  count 
retea  (raw)  In  IMF  positive  and  negative  aeetors 
aa  a  function  of  geoaagnetle  disturbance  (ap) 
for  the  years  196*  and  1965  of  eunspot  slniaua. 


Figure  2.  Mt.  Wellington  NM  data  (raw);  the 
aaae  conditions  described  in  Figure  1  apply  to 
all  *  figures.  The  high  noise  content  results 
primarily  fro*  the  low  count  rate  of  this  HM. 
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Figure  3.  Deep  River  DM  (filtered  to  attenuate 
the  statistical  noise).  Rotlce  the  value  of 
algal  given  Is  determined  by  counting  statistics 
only. 


Figure  ».  Mt.  Wellington  RM  (filtered).  In 
coeparlson  with  the  Deep  River  plots,  the 
saaller  eeparatlon  of  the  two  curves  and  the 
sealler  value  of  eexleua  slope  are  both  ascribed 
to  the  higher  cutoff  rigidity  at  Ht.  Wellington 
(see  text). 


Figure  Captions 


Figure  1.  Deep  River  neutron  monitor  (NM)  count  rates  (raw)  in  IMF  positive 
and  negative  sectors  as  a  function  of  geomagnetic  disturbance  (ap)  for  the 
years  1 96^4  and  1965  of  sunspot  minimum. 

Figure  2.  Mt.  Wellington  NM  data  (raw);  the  same  conditions  described  in 
Figure  1  apply  to  all  4  figures.  The  high  noise  content  results  primarily 
from  the  low  count  rate  of  this  NM. 

Figure  3.  Deep  River  NM  (filtered  to  attenuate  the  statistical  noise). 
Notice  the  value  of  sigma  given  is  determined  by  counting  statistics  only. 

Figure  4.  Mt.  Wellington  NM  (filtered).  In  comparison  with  the  Deep  River 
plots,  the  smaller  separation  of  the  two  curves  and  the  smaller  value  of  max¬ 
imum  slope  are  both  ascribed  to  the  higher  cutoff  rigidity  at  Mt.  Wellington 
(see  text). 


